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Abstract 

We consider a holographic description of the inflationary multiverse, according to which the wave 
function of the universe is interpreted as the generating functional for a lower dimensional Euclidean 
theory. We analyze a simple model where transitions between inflationary vacua occur through 
bubble nucleation, and the inflating part of spacetime consists of de Sitter regions separated by 
thin bubble walls. In this model, we present some evidence that the dual theory is conformally 
invariant in the UV. 
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I. INTRODUCTION 



Defining the probability measure in an eternally inflating universe is one of the key 
unresolved problems of inflationary cosmology. Eternal inflation produces an infinite number 
of "pocket universes" , in which all possible events happen - not once, but an infinite number 
of times. We have to learn how to regulate and compare these infinities, since otherwise we 
cannot distinguish between probable and highly improbable events, and thus cannot make 
any predictions at all. 

In "multiverse" models with a multitude of different vacua, eternal inflation gives rise 
to a fractal pattern, where pockets of all possible vacua are nested within one another. An 
important problem in this type of model is to find the probability distribution for the values 
of low-energy constants of nature, such as the parameters of the standard model of particle 
physics. Once again, we need some way of regulating infinities, and the problem here is 
exacerbated by the complicated spacetime structure of the multiverse. 

The essence of the problem is that the numbers of all kinds of events in an eternally 
inflating universe are growing exponentially with time. Whatever cutoff method is used, 
most of the events occur near the cutoff, and the resulting probability measure depends 
sensitively on the cutoff prescription. This is the so-called "measure problem" of inflationary 
cosmology. (For a review see l2|.) 

So far, most of the work on the measure problem has been phenomenological. Different 
measure proposals have been examined to check whether or not they lead to inconsistencies or 
to a glaring conflict with the data. (For recent discussion and references, see, e.g., 

mm 

.) Other selection criteria have also been introduced. For example, it has been argued that 
the probability measure should not depend on the initial conditions at the beginning of 
inflation. It seems rather unlikely, however, that this kind of analysis will lead to a unique 
prescription for the measure. 

A more satisfactory approach would be to motivate the choice of measure from some 
fundamental theor y. Attempts in this direction have been made in In 
particular, in Ref. [lO|] we suggested that the dynamics of the inflationary multiverse could 
have a dual description in the form of a lower- dimensional Euclidean field theory defined 
at the future infinity. The measure of the multiverse can then be defined by imposing a 
Wilsonian ultraviolet cutoff ^ in that theory. We argued that in the limit of 0, the 
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boundary theory becomes conformally invariant, approaching a UV fixed point. We also 
argued that on super-horizon scales the UV cutoff ^ corresponds to a scale factor cutoff in 
the bulk theory. 



In the present paper we shall further explore the holographic duality proposed in Ref . [10| , 
filling in some of the missing details. In Section II we consider the simple model where tran- 
sitions between different vacua occur through bubble nucleation, and the inflating part of 
spacetime consists of de Sitter regions separated by thin bubble walls. We discuss the struc- 
ture of future infinity in this model, focusing in particular on the eternal set defined by 
eternal timelike curves which always remain within the infiating region and never encounter 
terminal bubbles of negative or zero vacuum energy density. We show that (i) the boundary 
metric on £ can be chosen to be fiat and (ii) the bubble distribution on the boundary is then 
approximately invariant under the Euclidean conformal group, with the invariance becoming 
exact in the UV limit. This supports the conjecture that the dual boundary theory should 
be conformally invariant in the UV. 

By analogy with AdS/CFT correspondence 1^, Q, Q, we suggested in Ref. [l^ that 
the correspondence between the multiverse theory in the bulk and its dual on the future 
boundary is expressed by the relation 

^[0(x)] = j D(t) e'^l-^l = e'^t^l. (1) 

Here, S is the bulk action and the integral is over bulk fields (f) approaching the prescribed 
= 0(x) at the boundary. The amplitude \I'[0(x)] has the meaning of the wave function of 
the universe, and W[(f)\ is the effective action for the boundary theory with the appropriate 
couplings to the external sources 0. In Section III we describe this proposal in more detail, 
with emphasis on the IR/UV connection and the implications for the measure problem. 

To gain further insight into the properties of the boundary theory, in Section IV we cal- 
culate W for linearized perturbations around the model of nested bubbles. First we consider 
tensor modes in de Sitter space. This can be used to find correlators like {h{x)h{x.')) when 
the points x,x' are within the same bubble. The functional form of the corresponding W is 
consistent with that expected in a conformal field theory. Then, we consider fiuctuations of 
the bubble walls in the approximation where the self-gravity of the bubble can be neglected. 
In the boundary theory, the bubble walls mark the boundaries between regions with different 
central charge, or different number of field degrees of freedom. Such boundaries give a con- 
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tribution to the trace anomaly which depends on their shape. We evaluate this contribution 
from the asymptotic form of the bulk wave function for the bubble wall fluctuations. Once 
again, the result is consistent with conformal invariance of the boundary theory. 

Finally, in Section V we summarize our conclusions and discuss some open issues. 

When this paper was nearly completed, we learned of the work in progress by Stephen 
Shenker, Douglas Stanford and Leonard Susskind, which has some overlap with the ideas 
presented here. 



A related interesting development is the work of Freivogel and Kleban IJ] . They consider 



a model where bubbles nucleate in a de Sitter background, and compute correlators for 
operators which characterize the bubble distribution at the future boundary. They find that 
these correlators are conformally invariant, and they also discuss a dual CFT interpretation 
of their results. 



II. THE MODEL OF NESTED DE SITTER BUBBLES 



A. Flat foliations 



A diagram illustrating the causal structure of an eternally inflating spacetime is shown 
in Fig. [H Bubbles of all possible types nucleate and expand, rapidly approaching the speed 
of light. The worldsheets of the bubble walls can therefore be approximated as light cones 
in the diagram. The future boundary of this spacetime includes the singular boundary 
corresponding to the big crunch singularities of the negative-energy anti-de Sitter (AdS) 
bubbles, "hats" corresponding to the future null and timelike infinities of the Minkowski 
bubbles, and the eternal set £ which is the focus of our interest in this paper. AdS and 
Minkowski bubbles are called "terminal bubbles", since inflation completely terminates in 
their interiors. The eternal set £ consists of the spacelike future boundaries of the inflating 
de Sitter (dS) bubbles - or rather what remains of these boundaries after we remove the 
regions eaten up by terminal bubbles. We can think of £ as the set of "endpoints" of eternal 
timelike curves, which never encounter terminal bubbles. (A more formal definition will be 
given in Section III.) 

We shall start by considering a simplified model in which the inflating part of spacetime 
consists of pure dS regions separated by thin bubble walls. This part of spacetime can be 
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Minkowski 




bubbles. 



FIG. 1: Causal diagram of the inflationary multiverse. 

foliated by spacelike surfaces Sf, labeled by some coordinate t. The key observation now 
is that in our simple model of nested dS bubbles, the spacetime region of interest can be 



foliated by flat Euclidean surfaces. This is not difficult to understand IGj. A dS space 
inflating at a rate H can be represented as a hyperboloid 

X2 + r2 _ y2 ^ ^-2 (2) 

embedded in a 5D Minkowski space. The flat foliation of this space is obtained by slicing it 
with null hyperplanes 

Y -T = const. (3) 
The resulting metric can be written as. 

ds^ = -H-^dt" + e^'d^\ (4) 



5 



where t is the scale factor time, that is, time measured in units of H~^. (This coordinate 
system covers only half of the full dS space, but this is sufficient to cover the interior of any 
bubble.) 

The spacetime geometry for a bubble of vacuum V in a background of vacuum V can be 
obtained by matching the hyperboloid (jl]) representing the parent vacuum with a similar 
hyperboloid for the daughter vacuum, 

+ (r - YoY -T^ = H'-\ (5) 

where the displacement Iq depends on the tension of the bubble wall. The two hyperboloids 
are joined along a (2+l)-dimensional hyperbolic surface which represents the worldsheet of 
the wall (see Fig. [2]). 




FIG. 2: The spacetime representing a bubble of lower-energy vacuum (shaded) expanding in a 
higher-energy de Sitter space. One of the flat foliation surfaces is also shown. 
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Once again, a flat foliation can be obtained by slicing this hybrid construction with null 
hyperplanes ([3]). Each slice will then consist of a spherical region of vacuum V embedded 
in an inflnite region of vacuum V (or vice versa). Since the spatial geometry of the slices is 
flat both inside and outside the bubble, the volume that is removed from vacuum V is equal 
to the volume of vacuum V' which replaces it. As a result, the two regions match without 
any discontinuity, even though there is a (5-function curvature singularity along the domain 
wall in AD. 

The metric in the bubble interior can also be brought to the form (jlj), but since the 
expansion rates in the two vacua are different, the corresponding time variables do not 
match at the bubble wall. In a flat slicing, such as (jlj), the physical radius of the bubble is 
given by [igI 



where D = (if — /2q)^/^ and Rq is the constant intrinsic curvature radius of the worldsheet 
(this constant is determined by the bubble wall tension and the vacuum energies on both 
sides of the wall). The radius is the same when viewed from inside and from outside, and 
thus we must have 



For t,t' ^ 1, the time shifts in t and t' are related by Af = At and correspond to equal 
multiplicative changes in the scale factors a{t) = e*. 

This argument can be extended to any number of bubbles of any kind nucleating in 
the parent vacuum, as well as to bubbles nucleating within those bubbles, etc. A slice 
through the inflating region of the multiverse can then be chosen as a 3D surface having a 
flat Euclidean geometry. (The label t here can be chosen as the time variable in the parent 
vacuum.) Each bubble which is crossed by this surface is represented by a spherical region, 
with bubbles nucleating later having smaller images on . The image of an inflating bubble 
will generally have images of daughter bubbles nested within it. The image will therefore 
look like a "sponge", whose "pores" contain vacua of different types. Some of the "pores" 
may correspond to other inflating vacua, in which case they will themselves look like sponges, 
and so forth. The "pores" corresponding to terminal vacua are represented by "holes" in 
Sj. As we move to later and later times t, we will see more and more of this structure. In 
the limit of t ^ oo, each "sponge" becomes a fractal set, with almost all volume eaten up 



R^(t) = (H-^e' - D f + Rl, 



(6) 



t' = t + \n{H'/H) + 0{e-'). 



(7) 
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by terminal bubbles. 

We note that the fiat foliation we have just described is not the same as the scale factor 
foliation (with the usual definition of the scale factor through the expansion along a geodesic 
congruence). The reason is that geodesies undergo some focusing at domain walls and the 
congruence takes some time to adjust to the new expansion rate as it crosses from one 
vacuum to another. As a result, the fiat slices deviate from constant scale factor surfaces, 
typically by (5t ~ 1. This deviation, however, is rather insignificant, and we shall refer to 
the variable t as the scale factor time. 

The size distribution for bubbles of different kinds on the fiat slices can be found using 
the formalism developed in Refs. [l6|, ll7|. Each bubble can be labeled by two indices {ij}, 
where i refers to the vacuum inside the bubble and j to the parent vacuum in which it 
nucleated. The number of type-zj bubbles nucleated in a unit comoving volume, V{t) = e^*, 
in an infinitesimal time interval dt is given by 

dN^j = e''X,,f,H-'dt. (8) 

Here, Xij is the bubble nucleation rate per unit physical spacetime volume, fj{t) is the 
fraction of space in the 3D slice occupied by vacuum j, and there is no summation over j. 
The evolution equation for fi{t) can be written as 

dt 

where 



m 



and 



^ X^,-H-' (11) 



is the probability for a bubble to nucleate per Hubble volume per Hubble time. 
The asymptotic solution of at large t has the form 

f^(t) = f^'^ + s,e-^' + ... (12) 

Here, has nonzero components only in terminal vacua, —q < is the dominant eigen- 
value of the matrix M (that is, the largest nonzero eigenvalue) and sj is the corresponding 
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eigenvector.^ Substituting (fT2l) into ([8]), we have 

dNij = XySje^^-'^'H-^dt. (13) 

Now, consider a bubble of type ij which was formed at time t. The radius of this bubble 
at t' t is 

R{t',t)^Hj^e''-K (14) 
Expressing t in terms of R and substituting in f|T3|) . we obtain the bubble distribution 

dNij = CijE-^^-^^dR, (15) 

where Cij is a constant (-R-independent) coefficient. We thus see that the size distribution 
for all bubble types follows the same power law (fT5l) . 

B. Conformal invar iance 

It follows from Eq. f|T4|) that a forward time translation by At results in a rescaling of 
(large) bubble sizes by the same factor exp(At). Since the size distribution of bubbles (fT5|) 
is a power law, the form of the distribution is unchanged under the rescaling. 

Factoring out the scale factor e*', we can introduce the comoving bubble radius, 

r{t',t) = e^''R{t\t), (16) 

which approaches a constant value at t' — > oo, 

r ^ Hrh'\ (17) 

where t is the time of bubble nucleation. In terms of the comoving coordinates x, we thus 
have an asymptotically static distribution of bubbles at future infinity, 

dNij oc r-^^-i^dr. (18) 

This distribution was derived using the asymptotic solution f|T2|) for fj{t) at large t, so we 
expect it to apply only approximately, becoming exact in the limit r 0. The shape of 



^ It has been shown in jl7| that for an irreducible landscape of vacua, in which any vacuum can be reached 
from any other nonterminal vacuum by a sequence of transitions, the dominant eigenvalue is real, negative, 
and nondegenerate. 



the distribution at large r is not universal and is influenced by the initial conditions at the 
onset of inflation (e.g., the kind of parent vacuum we start with, etc.) 

The distribution (1181) is in agreement with scale invariance: in a space of dimension d, a 
scale invariant distribution has the form dN oc r^^'^^^^'dr, while the fractal dimension of the 




Another asymptotic symmetry of the bubble distribution is related to the possibility of 
choosing different flat foliations of spacetime. These foliations are related by Lorentz trans- 
formations in the 5D embedding space, or equivalently, by de Sitter group transformations 
in the parent vacuum region. Each foliation deflnes a coordinate system (jl]) and a con- 
gruence of timelike geodesies x = const in the parent vacuum. Any two such congruences 
become asymptotically comoving (in the region where they overlap) and thus deflne a coor- 
dinate transformation x ^ x' at the future inflnity. In the Appendix we show that it is a 
transformation of the form 

accompanied by a rotation. (Here, 6* is the dS boost parameter.) The transformation ( |T9|l 
is the so-called special conformal transformation (SCT). It can be described as a translation 
preceeded and followed by inversions. An important property of SCTs is that they map 
spheres into spheres, and thus the spherical shape of the bubbles is preserved. 

By the same argument, a change of foliation induces SCTs in the asymptotic future 
of all dS bubble interiors. By continuity accross bubble boundaries, it follows that all 
transformations should be the same, so that the entire future boundary £ is transformed by 
a single SCT (plus rotation). 

For any flat slicing, the evolution of the bubble distribution is described by the same 
equations (!8|) . (!9|) . (fT4l) . although the initial conditions are generally different. However, since 
the late-time behavior is universal and independent of the initial state, the form of the 
distribution should be the same at t — oo. This means that the size distribution of the 
bubbles should be invariant under SCTs in the limit of r ^ 0. 

Together with dilatations, translations, and rigid rotations, SCTs comprise the Euclidean 
conformal group. Our conclusion is thus that the bubble distribution is conformally invariant 
at r 0. If indeed the multiverse has a dual description on this suggests that the 
boundary theory should be conformally invariant in the UV. 
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C. Bubble collisions and the persistence of memory 



So far in this discussion we have ignored bubble colhsions. These are interesting in that 
they are sensitive to initial conditions, an effect which has been dubbed "the persistence 
of memory" 18|. This effect can be described as follows. Suppose we have some initial 
hypersurface of vacuum of type i which has no bubbles of any other phase. Then, memory of 
such initial condition persists arbitrarily far into the future. Indeed, the preferred congruence 
orthogonal to the initial "no bubble" hypersurface defines a rest frame. An observer in 
vacuum i at rest in this frame will be equally likely to be hit by a bubble of type j from any 
direction in space. In other words, the probability of collision with new bubbles is isotropic 
for this observer. However, if the observer moves with respect to the preferred congruence, 
then the observer is more likely to be hit head on by new bubbles. 

When we look at the shapes drawn by the nucleated bubbles at the future boundary, this 
translates into the following effect. Let us first choose our foliation to be (locally) orthogonal 
to the preferred congruence. In these co-moving coordinates, any bubble which nucleates 
in the original vacuum i will be equally likely to be hit from any direction. At the future 
boundary, this will be seen as a big spherical bubble "decorated" isotropically by a froth of 
smaller bubbles. 

However, if the foliation is boosted with respect to the preferred congruence, then the 
pattern of bubbles at the future boundary will be distorted by the corresponding SCT. 
Although the spherical shape of a bubble is preserved by inversions around an arbitrary 
point, the isotropic distribution of smaller spheres decorating it will be shifted preferentially 
in the direction determined by 6*. As a result, the froth decorating the bubbles will be 
anisotropic. 

Nonetheless, for a given value of the boost parameter 6*, the effect will get smaller for 
bubbles nucleating at later times. The reason is that a congruence orthogonal to any flat 
foliation becomes asymptotically comoving to the preferred congruence. In this way, the 
conformal invariance is recovered in the UV. 
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D. More general foliations 



Flat foliations of the kind we discussed so far are possible in the simple model of nested 
bubbles, but cannot be constructed in a general spacetime. To get an idea of what the 
general situation may be like, we shall now examine a wider class of surfaces. 

For a general spacelike surface S, we can introduce a coordinate system (jl]) in each 
dS region that this surface intersects. In terms of these coordinates, the surface can be 
represented as 

t = /(x). (20) 

The function /(x) may change its form from one bubble to another, but we shall require 
that S remains smooth at the bubble walls. We shall also assume that /(x) is a slowly 
varying function, so that 

e-^W|V/(x)|«i/(x). (21) 

This guarantees that the typical curvature radii of S are much greater than the local dS 
horizon. In other words, the surface is nearly flat on the horizon scale. 
The induced metric on E is 

ds^ = e^^^^^c/x^ - H-^dJdjfdx'dx^. (22) 

The condition ( l2Ti) allows us to drop the second term, so we can write 

ds^ ^ e^^Wrfx^ (23) 

We thus see that the metric on S is related to the fiat Euclidean metric by a Weyl rescaling 

H. 

In a CFT, the observables have well known transformation properties under Weyl rescal- 
ings. Once we know the observables for a given metric, we can find their values for metrics 
in the same conformal class. In this sense, the boundary theory is Weyl covariant, and the 
conformal factor is a gauge redundancy. We fix the gauge by choosing a particular metric 
in the conformal class. For instance, the choice of a fiat boundary metric has the advantage 
that the conformal invariance of the theory is manifest in some observables. 
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E. More general spacetimes 



To get an idea of what the situation is in more general spacetimes, not necessarily piece- 
wise de Sitter, let us consider the case of a metric which is locally FRW, but with the rate 
of expansion varying from place to place on scales much bigger than the horizon. To be 
more precise, let us consider the metric of a spacetime whose expansion is isotropic. This 
means that the extrinsic curvature of equal time slices orthogonal to a geodesic congruence 
is proportional to the spatial metric on the slices. The metric can then be written as 

ds^ = -dr^ + e^^jij{x)dx'dx\ (24) 

where N{t, x) can be interpreted as the number of e-foldings in proper time gauge. Intro- 
ducing r] = —e^^, we have 

dr] = -rjdN = -7]{HdT + N,i dx'), 

where H = N,t-. Hence, 

ds'^ = r]-^[~H-\dr] + r]N,i dx'f + ^ijdx'dx^]. (25) 

In the limit — 0, the second term within the round brackets can be neglected,^ and we 
have 

ds^ ^T]-^ [-H-^drf + -iijdx'dx^] . (26) 

Thus, we find that near the future boundary the metric of the inflating part of space-time 
is conformal to 

ds^ = -H-^{ri, x)dri'^ + -fij{x)dx'dx^ . (27) 

This suggests that we can identify S with the future conformal boundary of the metric (127!) . 
Clearly, jij plays the role of the metric at the conformal boundary. The arbitrariness in 
the choice of the congruence translates into Weyl rescalings of the metric at the confor- 
mal boundary. Indeed, in the inflating background, any two geodesic congruences become 
asymptotically co-moving to each other. Consider two geodesies that have the same end- 
point X at the conformal boundary. By the time their relative speed becomes completely 

^ In neglecting this term, we are implicitly assuming that the metric is smoothed over a fixed co-moving 
scale. The motivation for this smoothing will be clear from the discussion in the next Section, since the 
UV cut-off of the boundary theory corresponds to fixed co-moving scale. 
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negligible, the number of e-foldings from the initial surface for both geodesies will be differ- 
ent by some AA^(x), so in the asymptotic future we have i] = r]'e'^^^^^\ Hence, the spatial 
metric on surfaces of very small constant 7] will be related to that on surfaces of very small 
constant r]' by Weyl rescaling. 



III. THE WAVE FUNCTION OF THE UNIVERSE AND THE IR/UV CORRE- 
SPONDENCE 

The wave function of the universe can be expressed as a path integral 

s 

iS 



m= / D<Pe'^ (28) 

where the integration is over the bulk fields (f){x) interpolating between some initial condi- 
tions on hypersurface Sj and approaching the prescribed values 0(x) on the surface S, which 
plays the role of the future boundary of the spacetime region of integration. The nature of 
the initial conditions on Ej is the subject of some debate (for a review see {l9] and references 
therein), but it will not be important for our discussion here. The reason is that an eternally 
inflating region of spacetime quickly forgets its initial conditions. The memory of the initial 
state is retained only on the largest comoving scales, while on smaller scales the evolution 
exhibits an attractor behavior. For example, we saw in Section II. A that the asymptotic 
bubble distribution (fTSll at short distances is independent of the initial conditions. This 
distribution is reached in any comoving region including at least one eternal geodesic. 

Thus, to study the asymptotic small-scale behavior, we do not need to invoke the wave 
function for the entire universe. We could start, say, with a cubic region filled with vacuum 
i and having size of a few H~^. The conditions at the boundaries of this region are not 
important; we could, for example, impose periodic boundary conditions. In our simple 
model, the spacetime consists of nested regions of dS space separated by thin walls, and the 
evolution is nearly classical, except for occasional nucleation of bubbles. According to the 
discussion in Section II. A, this spacetime can be foliated by fiat hypersurfaces, so the initial 
and final surfaces can be chosen to be flat. The future boundary can then be approximated 
by one of the surfaces and the only variables that need to be specified on that boundary 
are the types of vacua and the centers and radii of the corresponding bubbles. We expect 
that in the limit t ^ oo all correlations of physical significance will not depend on the choice 
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of the initial vacuum on Ej. 

The conjectured correspondence relates the wave function for the bulk theory to 
the effective action of the boundary theory. The surface S, with regions covered by terminal 
bubbles removed, gives an approximate representation of the eternal set S, as we describe 
below. As we shall see, the boundary theory, regularized with an appropriate UV cutoff, 
can also be thought of as living on that surface. 

A. The discrete boundary theory and its continuum Hmit 

A point P E S can be thought of as the endpoint of an eternal timelike curve (that is, 
a curve that never encounters terminal bubbles). More precisely, "points" or elements of £ 
are identified with chronological pasts of eternal timelike curves [20]. Two curves with the 
same past define the same element of S. Physically, an element of S is the equivalence class 
of eternal time-like curves which remain forever in mutual causal contact. An eternal curve 
can be pictured as the worldline of an observer and the corresponding boundary "point" P 
is the spacetime region inside the causal horizon of that observer (see Fig. [3]).^ 

The image Ps of a point P E S on the surface S is given by the intersection of S with P 
(see Fig. [3]). In a pure dS space, Ps would be a spherical region of radius equal to the dS 
horizon. In our model, the shape of Ps will be distorted by intervening bubbles, but its size 
will still be roughly given by the local horizon. (If Ps covers more than one vacuum, its size 
is typically set by the largest of the corresponding horizons.) In this way the entire eternal 
set S is mapped onto S, with each point represented by a roughly horizon-size region. 

Now, instead of considering the boundary theory at S, let us consider a boundary theory 
defined on S, which should approximate the theory at S on scales larger than a certain 
Wilsonian cut-off. 

First of all, we note that from the point of view of the theory living on E, the maximum 
possible resolution is the size of the images Pe of the points P of the future boundary. This 

Elements of £ can also be identified with points in the spacelike part of the future conformal infinity, like 
the point P in the conformal diagram in Fig. [3l This point is the endpoint of the eternal timelike curve 
7, and the past of 7 appears in the diagram as the interior of the past light cone of P. This definition of 
£ is possible only if the spacetime admits a conformal infinity, which may not be the case in general. 
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FIG. 3: Eternal curves 7 and 7' have the same past and thus define the same point P £ £. The 
intersection of this past (appearing in the figure as the past hght cone of P) with a spacelike surface 
S gives the image of P on S. 

resolution is position dependent, 

/,(x) ~ H'\x), (29) 

where H~^{x) is the local dS horizon radius. It is interesting to note that, with this choice, 
the number of degrees of freedom of the boundary theory on S corresponds to that of the 
3ulk theory."^ The argument here is similar to that of Susskind and Witten for AdS/CFT 

1- 

The number of bulk degrees of freedom in a horizon region of vacuum i is given by the 



2l|. 



^ In the present discussion, and following standard practice, number of degrees of freedom is synonymous 
with the logarithm of the number of quantum states that the system can be in. 
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Gibbons-Hawking entropy of the dS space, 

M^uik ~ H-\ (30) 

To estimate the corresponding number Mboundary in the boundary theory, we note that, with 
the cutoff scale aX Ic = 1/H, the energy density of each field at the boundary cannot much 
exceed 1 / H^. This means that a horizon volume can contain at most ~ 1 quanta of frequency 
~ H. Longer wavelength modes contribute much less to the entropy, so we can disregard 
them. Thus, each field carries about one bit of information per horizon, and the entropy per 
horizon is of the order of the number of fields.^ For a conformal field theory, this number 
is related to the central charge of the theory. For the CFT dual to dS space, it has been 



estimated as 



22 



23| 



■^boundary ~ ^. (31) 

(see also the discussion in Section IV.) Thus, we see that, with Ic ~ H^^, we have 

Mbulk ^ -^boundary (32) 

The theory on E can be thought of as a field theory only on scales larger than /c, which 
plays a role analogous to the lattice spacing in a discrete system. Borrowing the terminology 
of Ref . , we can call this irregular lattice "a fish net" . Each vertex of this fish net carries 
■^boundary "internal" degrees of freedom. 

The continuum limit on S is obtained for physical wavelengths longer than a certain 
Wilsonian cutoff ^o, which should be much bigger than all "lattice spacings" H~^^ 

^0 > ^0 = maxiijr\ (33) 

where the maximization is over all dS vacua. In the nested dS model, we can choose E to 
be one of the surfaces S^; then the Wilsonian cutoff ^ in co-moving coordinates must satisfy 

i > L{t) = Lq e-\ (34) 

Here, L{t) can be thought of as a physical UV cutoff, where the continuum limit breaks 
down. This cutoff shrinks as the bulk boundary surface S is moved to the future. As a 
result, at later times the eternal fractal 8 can be seen at a greater and greater resolution. 



We thank Raphael Bousso for clarifying this point to us. 
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This situation is analogous to the IR/UV correspondence in the AdS/CFT [2l|, with the late- 
time boundary playing the role of the IR cutoff in the bulk. The scale-factor time evolution 
in the bulk corresponds to the renormalization group flow in the boundary theory.^ The 
theory on £ is obtained in the limit t ^ oo, where ^ and there is no physical UV 
cutoff. 

In our discussion so far we assumed the spacetime to the future of S to be well defined 
and classical. Of course, this is not so, even in the nested dS model. Spontaneous nucleation 
of bubbles allows for a multitude of possible histories, and the mapping of £ onto S will 
generally be different for different histories. We note, however, that these differences affect 
the mapping only on scales smaller than or comparable to the causal horizon, which in turn 
is smaller than Lq. So, with the choice of the cutoff length satisfying fl5^ . the mapping is 
insensitive to the future evolution. 

The picture suggested by the above discussion is the following. The hypersurface E can 
be thought of as a discrete system consisting of a juxtaposition of "lumps" Ps- From a 
bulk point of view, a lump is the interior of the causal horizon of an eternal observer. This 
has a finite number of internal degrees of freedom, proportional to the horizon area (it may 
therefore be appropriate to visualize as a horizon sized closed shell embedded in S, with 
the degrees of freedom living on the shell). At distance scales much larger than the size of 
the lumps, a continuum description emerges. The co-moving size of the lumps gets smaller 
as the reference surface S is pushed forward in time, and in the limit when it is sent to 
future infinity, the continuum description becomes valid all the way to the UV (where it 
should be conformally invariant). 



B. The boundary measure 



In Ref. [lOl] we proposed that a measure for calculating probabilities in the multiverse 
may be formulated at future infinity £. The idea was to use the Wilsonian cutoff ^ of the 
boundary theory to render the number of events finite. Any event occurring in the bulk, 



6 



For the case of Euclidean AdS, the CaUan-Symanzik RG flow equations for the boundary theory can be 



derived from a semiclassical expansion of the Wheeler-de Witt equation in the bulk [2j|, with the scale 
factor playing the role of the renormalization scale. Similar considerations were presented in Refs. 25, 2^ 
for the case of slow-roll inflation. 
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which requires for its description a bulk resolution corresponding to a physical wavelength 
Amin, will leave an imprint on £. The resolution which is needed to reconstruct the event 
from its imprint at the future boundary will be given by 

C — ^min 6 ** , (35) 

where is the scale factor time of the event. For a fixed bulk resolution Xmin, the Wilsonian 
UV cutoff ^ of the boundary theory corresponds to an IR scale factor cutoff at t = t* in the 
bulk, with ^ and t^: related by Eq. ( l35l) . Note that Xmin, which plays the role of a Wilsonian 
UV cut-off in the bulk, depends in principle on what kind of event we are interested in, 
and the cutoff time will have that dependence as well. An alternative approach is to set 
Xmin ~ 1, so that all events above the Planck scale will be resolved at t < t*. 

Our prescription for the boundary measure is not limited to the nested dS model or flat 
foliations. In the general case, we can use a scale factor foliation, starting with some surface 
So, which is smooth on scales larger than Lq and otherwise arbitrary. The scale factor a is 
defined as the expansion factor along the congruence of geodesies orthogonal to Sq, and the 
scale factor time is defined as t = In a. As before, each surface of constant t is the site of 
a boundary theory with a UV cutoff ^ satisfying flMl) . and the scale-factor time evolution 
corresponds to RG flow in the boundary theory. 

The scale factor time is not a good foliation parameter in regions of structure formation, 
where geodesies converge and cross. However, these phenomena affect only sub-horizon 
scales, which are smoothed out by the super-horizon cutoff ^. We require, therefore, that 
the foliation surfaces should also be smooth on scales ^ ^. 

Different choices of the initial surface Sq are related by Weyl transformations on S. The 
freedom of choosing this surface can be used to obtain a foliation with desired properties. 
A flat foliation, like the one we discussed for the nested dS model, is possible only in very 
special cases. Its closest analogue in a more general spacetime is a foliation by surfaces 

^ The correspondence between the boundary cut-ofF and the scale factor cut-off is only approximate. The 
main reason is that the scale factor time is not well defined on subhorizon scales. If one uses the density of 
a dust of test particles in order to define the expansion, then this becomes ill defined in regions of structure 
formation. One can use other definitions, but the choice is not unique. Moreover, Eq. ([35)1 assumes that 
wavelengths are conformally stretched with the expansion. This is certainly the case on superhorizon 
scales, but on subhorizon scales the wavelengths of signals can be affected by causal processes other than 
the expansion of the universe. 
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having a vanishing Ricci scalar, 

R^^^ = 0. (36) 

In a spacetime which is locally FLRW with small perturbations, this condition can always be 
satisfied by a position-dependent time shift. If we choose the initial 3-surface Eg satisfying 
( l36l) . it can be shown that this property is preserved (to linear order) by scale factor time 



evolution 
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IV. THE WAVE FUNCTION AND CONFORMAL INVARIANCE 

In this Section we go beyond the simple model of nested bubbles discussed in Sec. II and 
allow some perturbations in the infiating regions and on the spherical bubble walls. First, 
we consider massless fields (or metric perturbations) in the case of a single vacuum. Then, 
we consider the bubble walls which separate domains with different vacua. In general, the 
bubble wall fluctuations are entangled with the metric perturbations. Here, we shall consider 
bubble fluctuations in the region of parameter space where the self-gravity of the bubble 
can be ignored and the two types of perturbations can be disentangled. As we shall see, 
the effective action which is obtained through the correspondence ([1]) does have some of the 
features expected in a CFT. 



Fort 
in Ref . 



le case of massless flelds in de Sitter, our discussion closely follows that of Maldacena 



28| . A formal diference is in the expression of the conjectured dS/CFT correspon- 
dence. We have on the right-hand side of ([1]), while Ref. 28| has . If the effective 
action is expressed in terms of a boundary path integral, we propose that this should be 
written in the form 

e^m^^ = j^D^e'^^^'^\ (37) 

There is no time variable at the boundary, and in this sense the action S for the boundary 
flelds ip is Euclidean. However, instead of using , as in standard Euclidean theories which 
are obtained by Wick rotation from Lorentzian time, here we propose using the phases e*"^. 
This is because the wave function \E^, given by the bulk path integral in ([1]), is complex.^ 
In fact, the phase of the wave function grows at late times in proportion to powers of the 
scale factor. In the boundary theory, these growing phases can be interpreted as ultraviolet 



This is in contrast with the Hartle-Hawking wave function, which is real. 
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divergences, which can be removed by local counterterms in 5*. With our conventions, the 
counterterms are real. 



A. Massless fields in de Sitter 

Let us consider the wave function of a free massless scalar field h in de Sitter. The same 
wave function describes linearized gravitons in the transverse traceless gauge. For the sake of 
comparison with existing CFT calculations, let us work with a (d+l)-dimensional de Sitter 
space. In flat conformal coordinates, the metric is given by 

ds^ = (ri) [-drf + c/x^] , (38) 

with a{ri) = —1/Hr]. Decomposing the field in Fourier components 



e 



ikx 



Mx) = J d'kj^^h^, (39) 
the Gaussian solution of the Schrodinger equation takes the form \E' = e*^, with 

jdu f^'^ ^ '"kii, |2 



W = J d'^k ^|/ikr + ^lnWkJ • (40) 

Here a prime indicates derivative with respect to r], and Vk{ri) is a solution of the massless 
wave equation v'^+ {d — l){a'/a)v{^ + k^fk = 0, which should be of "negative frequency", 
^k^k ~ ^k^k — so that the Gaussian is normalizable. The Bunch-Davies vacuum 

corresponds to the choice 

MV) = ^,^-"\v)H^lUkr^). (41) 



where Hankel function. 



Let us consider the case d + 1 = 5. At late times, after modes have crossed the horizon 
{—krj <^ 1), Eq. (HOl) takes the asymptotic form 

WU^)] = \ |rf^k(^ + ^[ln(A:V/7V)+m + 27]+0(a-2))|/.kP + ..., (42) 

where the ellipsis denote terms which are independent of h\^. The structure of (H2l) is similar 
to the one obtained in the context of AdS/CFT by Gubser, Klebanov and Polyakov fl^ . 
except that, here, W contains an imaginary part. The imaginary part should be there 
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because the amplitude of perturbations on superhorizon scales is determined by 
g-2/m[VK]_ 'pj-^jg ]-^g^g y^oii defined limit as a oo, 



oc exp 



71 



(43) 



corresponding to a scale-invariant spectrum of perturbations {hl.h-\^i) = {8H^ /TTk^)6(k.' — k). 
Let us now consider the real part of W. The first term diverges like as we approach the 
future boundary. This term is analytic in k"^, and its coefficient can be changed by adding 
boundary counterterms of the form J {dihYd^^. The divergence of the kinetic term as a — cx3 
suggests that the field is non-dynamical at the boundary. Introducing /i = — r/^^ = Ha, the 
non-local part is written as 

Re[W] = ^ [ d^k fc^ln(A;V/i2) \hk\^ + analytic. (44) 
16 J 

As mentioned above, the wave function for linearized gravitons in the traceless and trans- 
verse gauge is exactly the same as that for the massless scalar field h. The expression f lU]) 
we obtain for W does indeed take the standard form of an effective action for a conformal 
theory coupled to an external gravitational field in 4 dimensions. For the case of free fields, 



this was first discussed by Tomboulis 
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The structure of (jHj) is not difficult to understand. The renormalized effective action for 
a CFT propagating in a curved background takes the form [s^ ^ 

Wren = ad/2 lu /i^ + . . . (45) 

where the ellipsis indicate terms which are independent of /x, and 

ad/2= [ {T)^d''x (46) 



is the integrated trace anomaly in d dimensions. When d is even, ad/2 is given in terms of 
an integral of geometric invariants constructed from contractions of the Riemann tensor. In 
d = 4, this takes the form 

a2 = J d^x^ [ciR^ + C2RijR'' + c^RijkiR''"'] • (47) 



^ This can be understood as follows. A change in the renormalization scale fi ^ fi is equivalent to a global 
rescaling of the metric gtj ^^gij, with ft = fl/fj,. Hence dWren/dlnfi = / {SWr en /S In n)d'^x. Since the 
trace anomaly (T) = {1 / y/ — g)SWren / S In should be independent of /i, it follows that Wren is linear in 
log/Lt, which leads to the form (|i5)) . We thank Igor Klebanov for a discussion of this point. 
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The coefficients q are such that 02 is a hnear combination of the Euler number and the 
integral of the Weyl tensor squared. Both are invariant under Weyl rescahngs (see e.g. 
{30! and references therein). More generally, 0^/2 is Weyl invariant in d dimensions. ^'^ To 
quadratic order in the metric perturbation hij, this leads to the structure Cj /c^j/ikP In/i^ in 
momentum space. The form of the non-analytic piece in is recovered by noting that 
the argument of the logarithm must be dimensionless. 

The numerical constants Cj in front of the geometric invariants scale in proportion to the 
number of fields (central charge) in the CFT. Comparing with (jUj) we have 

c ~ H-^. (48) 



Parametrically, this is also the entropy of the bulk de Sitter space [22|, l23|, l28 |. 

In the case of a 4 dimensional bulk {d + 1 = 4), substituting (HTj) into (HUjl . and expanding 
at late times {—kt] -C 1), one obtains 

= \ jd'k + + 0(a-^)) l/^^r + (49) 

The first term, which diverges as we approach the future boundary, is analytic, while the 
second one is non-local and imaginary. The second term leads to the familiar expression for 
the scale invariant amplitude of tensor modes (/ij^/ik') = H'^k~^6(k' — k). In the present case, 
the Hankel functions are of half integer order, and do not lead to logarithmic terms in the 
asymptotic future. From the CFT point of view, this is in agreement with the fact that there 
is no trace anomaly in odd dimensions, and so there are no logarithmic divergences. The 
number of fields in the CFT can be estimated from the coefficient of the non-analytic piece 



in W. Again, this scales in proportion to the entropy of the bulk de Sitter space [22|, l23l. |28| . 

c ~ H~^. (50) 



B. Bubble fluctuations. 



Throughout this subsection we work in the approximation where gravity of the bubble is 
unimportant. In the thin wall limit, we require that 

TRo < 1, iApv)Rl < 1. (51) 



This follows from (|45l) . by noting that the functional derivative of Wren with respect to Inil is the 
anomalous trace (T), which should not depend on /x. 
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where T is the tension of the bubble wall, Apv is the energy density gap between the inside 
and the outside of the bubble wall, and Rq is the intrinsic curvature radius of the worldsheet. 
The conditions flSlI) mean that the geometry on the scale of the horizon is not appreciably 
distorted due to gravity of the bubble wall and due to the change in the vacuum energy. 



Under these approximations, Rq is given by 3l|, |32| 



R^ ~ (^+^)^^^ (521 

For reference, we have written the expression for arbitrary worldsheet dimension p + 1, 
although we are mostly interested in the case of ordinary membranes p = 2. Note that we 
allow the energy gap Apv to be significant, even comparable to the energy density in the 
false vacuum, as long as ( ISTl) is satisfied. 

The interest of this weak backreaction limit is that we can study fiuctuations of the 
bubble wall without the need of coupling these to metric fiuctuations. We will call this 
the "Goldstone limit", since the only relevant degree of freedom is the one associated with 
normal displacements 6x^ of the worldsheet a;^(^*). In the fiat chart of de Sitter space, 

ds^ = -dt^ + e^H\dr^ + r''dn'') (53) 

the trajectory of a bubble wall centered at the origin of coordinates is given by 



(t) = /J-2(e-2^*« + - 2H-\H-^ - i?^) VSg-i/im,)^ (54) 



where to is a free parameter which can be interpreted as the time of nucleation. The unit 
vector normal to the wall worldsheet is given by 



r 

n 



i?o-' {{H-^ - Rlf'^ - i/-^e-^(*-*°)) , (55) 
Ro^r^{t)e-''^'^'''\ (56) 



We now parametrize the normal displacement of the wall as 33 1 



The factor T"^/^ jg introduced so that <^ is a worldsheet scalar with standard normalization. 
The linearized action for this scalar consists of a canonical kinetic term and a tachyonic 
mass term: 

m5 = -(p + l)i?o'- (58) 
24 



Let us now study the wave functional for this scalar field 0. The worldsheet of the bubble 
is itself a de Sitter space of curvature radius Rq and dimension p + 1. This worldsheet will 
cut a sphere at future infinity, and for that reason it is convenient to use the closed chart for 
the p + 1 dimensional de Sitter spacetime, with coordinates ^* = {fj,Qa), with a = 1, ...,p, 
and metric given by 

dsl = {-df + dnl) . (59) 

Here 

d = Rq/ cosfj, (60) 
and dQp is the metric on a unit p- sphere. Expanding the normal displacement as 



J2^LMYLM{n), (61) 



LM 

where Yla/ are spherical harmonics, the wave functional is given by ^ = e*^, with 

W = J2 S<PLM? + ^\nvI\ . (62) 

The Bunch-Davies vacuum in the closed chart corresponds to the choice 



vl 



ALicosvy/' (^P^i+^/2(sin77) + ^Ql_,^^/,ismv)^ , (63) 



where is a normalization constant. The index of the Legendre functions P and Q is given 
in terms of the mass of the field as 

z/ = I — - m^i?o 1 = (64) 

For bubbles in the 3+1 dimensional multiverse we have p = 2, and u = 2. Inserting (163|) 
into (|62|) and expanding for late times (a oo, or fj —>■ 7r/2), we have 



LM 



2Rn ' 4 ' 16 52 



ln(||)+2(^(L) + l)+.vr + 27 



-'LA/ 



|2 



(65) 

where we have omitted terms of order (as well as terms independent of (p). Here, is the 
digamma function, and we are using A = —L{L + 1) to denote the Laplacian eigenvalues. 



To simplify this expression, we have used a momentum independent field redefinition at late times (p = 
{sin fiY^'^cf), when sinry = 1 + 0(a~^). In the expression above we have not written the hat on top of (f>LM, 
since we will be interested in the behaviour at the future boundary, where both variables coincide. 
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In the limit of large L, ip{L) ^ InL, so the non-analytic term in (l65i) is proportional to 
ln(i?QL^/4a^). The formal similarity between fH2l) and fl65l) is to be expected, since the 
index z/ = 2 of the Legendre functions which characterizes a field with m?R^ = — 3 in the 
2+1 dimensional de Sitter worldsheet also corresponds to a massless field in 4+1 dimensional 
de Sitter space, considered in the previous subection. For that reason, the mode solutions 
are very similar (except for an overall power of the scale factor). 

Note that the imaginary part of W , contributing to l^&P, decays as a~^. This is not 
surprising, since the normal displacement is a tachyonic field whose amplitude grows with 
the scale factor (pLAiiv) oc a(f/) at late times. 

To make contact with the boundary theory, we would like to consider the relative co- 
moving displacement on the fiat slices of constant scale factor time. This is given by 

The second factor in the right hand side is the proper normal displacement (in the reference 
frame where the worldsheet is at rest) divided by the physical radius of the bubble. The 
relativistic 7 factor accounts for the Lorentz contraction of the displacement in the reference 
frame associated to the constant scale factor hypersurfaces, and can be calculated from fl56|) 

as 

J = n^^n^ = R-'r^{t)e"'^ {HR^)-\ (67) 

Here, = a{t)5^r is the normal to a co-moving sphere, and in the last step we have taken 
the late time limit. 

Note that 7 will be different on both sides of the wall if the vacuum energy is different 
on both sides. Geometrically, this means that the surfaces of constant scale factor time 
on both sides of the wall meet at an angle, and the corresponding frames are not at rest 
with respect to each other. Although the unperturbed geometry of the surfaces is 
fiat, the embedding of these surfaces in the bulk is non-trivial. Consequently, the relative 
perturbations measured by co-moving observers on both sides are related by 

= H;'6,, (68) 

Here we are neglecting the gravitational field of the domain wall, so this kinematic effect is unrelated to 
the jump in the extrinsic curvature as we go accross the worldsheet, which would only be appreciable if 
the inequalities ([5T|) were violated. 
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where the indices i and e refer to the interior and the exterior. This means that in the 
presence of bubble fluctuations, we cannot represent the future boundary as a smooth flat 
surface, because the interiors of the perturbed bubbles do not fit nicely within the holes 
they carve in the parent vacuum (the relative size of the wiggles differs on both sides). On 
the other hand, the bulk geometry is smooth, and in principle we could choose to foliate 
the perturbed bubble smoothly. In this case, the interior of the bubble will match the 
exterior, even as we approach the future boundary, but the metric on the surfaces S' of the 
new foliation would not be flat: there would be some metric perturbations. In this sense, 
the situation would not be so different from the case considered in the previous subsection, 
where we allowed for metric perturbations on a de Sitter background. 

For the remainder of this Section, we shall avoid introducing curved foliations by restrict- 
ing attention to the case where the jump in energy density accross the wall is sufficiently 
small, 

AH , , 

« 1- (69) 

In this limit, we can think of the bubble as propagating on a fixed background de Sitter 
space characterized by a single parameter if, in which case 5e ~ ^i- Using (l66l) and (1671) . 
Eq. ([HS]) can be written as 



5la/|' + ... (70) 



LM 

where we have dropped the analytic divergent terms. 

The expression vanishes both for L = (A = 0) and L = 1 (A = —2). This is in 
agreement with conformal invariance, which requires that the effective action is independent 
of rescalings of the bubble size and linearized translations (corresponding to L = and 
L = 1 respectively). 

As mentioned in the previous subsection, for a CFT propagating in a curved background, 
the coefficient of the logarithmic divergence is equal to the integrated trace anomaly co- 
efficient When d is odd and the manifold has no boundaries, a^j^ vanishes because 
we cannot build curvature invariants of odd dimension. However, if the manifold has sur- 
faces of co-dimension 1 where the CFT fields satisfy boundary conditions, then there is a 
contribution to a^i2 from geometric invariants constructed out of the intrinsic and extrinsic 
curvatures, Rahcd and Kah-, of the p-dimensional boundary (where p = — 1), as well as 
contractions of the (i- dimensional curvature Ri^ with the normal to the boundary surface. 
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If the boundary conditions are Weyl invariant, then the corresponding coefficient 0^/2 is also 
Weyl invariant. 

Let us now argue that the locus of the bubble wall at the future boundary plays the role 
of such a boundary surface. The case of our interest is p = 2, where the invariants have to 
be of dimension 2, and so they must be linear in Rijifn^ or R, or quadratic in the extrinsic 
curvature Kij. If Qij is fiat, as we are assuming in this subsection for the metric at the future 
boundary, then Rij = 0, and so the integrated trace anomaly must be of the form 



03/2 



rfi ( KabK''' - 1 + d2R 



(71) 



Here, dll2 is the area element on the boundary surface and di and d2 are constants. It can 
be shown that the term accompanying di is Weyl invariant, while the term accompanying 
d2 is topological. It is easy to check that in linearized theory around a spherical bubble wall, 
we have 

j dT.2 (^KabK"^ - ^K^^ (X J dn 6A{A + 2)6, (72) 

where 6 = 6r/ru, is the relative radial displacement at the future boundary, and f2 is the 
solid angle. Note that the coefficient of the logarithmic divergence in fITOl) has precisely this 
form. On the other hand, since the term proportional to d2 is topological, it does not depend 
on the perturbation 6. 

To conclude, the logarithmic divergence in fl70p takes the form expected in a CFT at the 
future boundary, where the domain walls play the role of surfaces where the CFT fields have 
to satisfy boundary conditions. 

The idea that domain walls act as boundary surfaces seems very natural in the present 
context, since the number of field degrees of freedom c ~ is different on both sides of 
the wall. The coefficient di can be read from (ITOil . 

~ ^- (73) 

Unless there are drastic cancellations, di should be roughly equal to the number of fields 
satisfying boundary conditions on the wall. It follows from f lSTl) that di ^ H^^, indicating 
that most of the fields pass freely through the bubble wall. 

On the other hand, the number of fields satisfying nontrivial boundary conditions on the 
wall should at least be equal to the difference Ac in the number of fields on the two sides of 
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the wall. In the absence of dramatic cancellations, this leads to 

c/i > Ac ~ AH/H^, (74) 
Substituting (!73l) into fl74l) . we obtain 

TRo > AH/H. (75) 

Using (1321) . the above inequality can only be satisfied provided that 

Apv < HT. (76) 

For such values of the parameters, the intrinsic curvature radius (or inverse proper acceler- 
ation) of the bubble walls is comparable to the inverse Hubble radius of the parent vacuum 

Ro ~ H-\ (77) 

The necessity of the condition (!76|) . and its possible implications, are at present unclear 
to us. It could be that there are unexpected cancellations in the boundary theory which 
allow di ^ Ac. On the other hand, it should also be noted that when (!76l) is violated, 
the proper acceleration Rq^ of the bubble walls is much larger than H, all the way to the 
future boundary. It could be that the correspondence between bulk and boundary theory in 
this "high energy" regime is not as straightforward as it seems to be for low energy bubbles 
(with Rq^ ~ H). The investigation of this issue is left for future research. 



V. CONCLUSIONS AND DISCUSSION 

We have explored the conjecture, made in Ref. [l^, that the infiationary multiverse has 
a dual holographic description at its future boundary, in the form of a lower dimensional 
theory which is conformally invariant in the UV. The duality is expressed by the relation 

^ = e*^, (78) 

where is the wave function of the multiverse with arguments on a spacelike hypersurface 
S and W is the effective action of the boundary theory on E. Here, we have argued that 

Supersymmetry will not necessarily help enforcing the cancellations which would make di much smaller 
than the number of fields satisfying boundary conditions at the wall. For instance, in N=4 SYM in d = 4 
we have supersymmetry, but the coefficient of the trace anomaly is comparable to the number of fields. 
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the boundary theory is defined on a "fish net" with characteristic spacing set by the local 
horizon size, and that the number of degrees of freedom in this theory is comparable to that 
in the bulk theory. A continuum description is obtained by imposing a Wilsonian cutoff, in 
the limit where the cutoff length C, is large compared to all fish net spacings. 

To study the UV limit of the boundary theory, we foliate the bulk spacetime (excluding 
terminal bubble interiors) with surfaces of constant scale factor time t, starting with some 
initial surface Sq. (Geodesic crossings that may occur in structure formation regions on 
sub-horizon scales do not interfere with this construction, since the foliation surfaces are 
smoothed out with a super-horizon cutoff ^.) As we go to larger values of t, the comoving 
fish net scale decreases as e~*, so we can choose the cutoff ^(t) oc e~*. Thus, renormalization 
group flow in the boundary theory corresponds to scale factor time evolution in the bulk, 
with the UV limit ,^ — > on the boundary corresponding to the IR limit t ^ oo in the bulk. 
In this limit, the foliation surfaces approach the eternal set S at the future inflnity. Different 
choices of the initial surface Sq are related by Weyl rescalings in the boundary theory. They 
should be physically equivalent in the UV if the theory is indeed conformally invariant in 
that limit. 

To find evidence for (or against) the conformal invariance of the boundary theory, we 
have studied a simple model in which the infiating bubble interiors are pure de Sitter, so 
the infiating part of spacetime consists of de Sitter regions separated by thin bubble walls. 
In this case the foliation surfaces can be chosen to be fiat, and we found that the bubble 
distribution on these surfaces is approximately invariant under the Euclidean conformal 
group, with the invariance becoming exact in the limit t oo. 

We have also studied the effect of linearized perturbations about the model of nested dS 
bubbles. Using the duality relation fl78l) . we have calculated the effective action W for the 
case of linearized tensor modes in de Sitter space and for fluctuations of the bubble walls 
in the limit in which the gravity of the wall is unimportant. In both cases, the resulting 
functional form of W is consistent with that expected in a conformal fleld theory. We 
interpret the locus of bubble walls at the future boundary as deflning the surfaces where 
CFT flelds must satisfy boundary conditions. The form of the logarithmically divergent 
terms in W is in agreement with this interpretation. Altogether, our results support the 
conjecture that the boundary theory is conformally invariant in the UV. 

A puzzle arises for the case of bubble walls whose proper acceleration Rq^ is much larger 
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than the Hubble rate H of the parent vacuum. In this "high energy" regime, the numerical 
coefficient di in front of the logarithmic divergence in W is much smaller than the number 
of CFT fields which must satisfy boundary conditions at the wall. This result is somewhat 
counterintuitive, and its implications are left for further research. 

Another important issue that has not been addressed in this paper is the treatment of 
terminal bubbles. We simply adopted the proposal of Ref. [lo| that the interiors of such 
bubbles should be excised on the future boundary, with the expectation that their dynamics 
is holographically described by 2d conformal field theories living on the bubble boundaries. 
Freivogel et. al. [7] provided some evidence for this in the case of Minkowski bubbles. For 
an AdS bubble, most of the volume in the interior is near the bubble wall, and it is natural 
to expect that a holographic description should apply. But as of now, this conjecture is not 
supported by any quantitative evidence. 

An objection that has often been raised against a holographic description of de Sitter 
space is the anomalous behavior of massive fields in the boundary theory. Strominger 22 1 
has studied the asymptotic behavior of the two-point function for a massive scalar field at 
future infinity and concluded that the field acquires a complex conformal weight if its mass 
is m > 3if/2. We note, however, that there may be a good physical reason for this behavior. 

Fields with m > 2iH/2 describe massive particles whose density is diluted as n oc a~^, so 
the comoving number of particles is conserved in the absence of interactions. The particles, 
however, are unstable and will decay in dS space, even if they are absolutely stable in the 
flat space limit: the effective Gibbons-Hawking temperature of the dS space allows decays 
in which the energy of the decay products is higher than that of the initial particle. Such 
particles cannot propagate any information to future infinity, so there is little to be learned 
from the asymptotic properties of their two-point function. On the other hand, the comoving 
number of particles with m < 3i7/2 grows with the scale factor time as 

oc e^* (79) 

with 



2m ^ ^ 



1 . 

3H 



1/2 



(80) 



If their decay rate is not too high, such particles have a chance of producing an imprint at 
future infinity. 
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As we mentioned in the Introduction, one of the main motivations for studying the 
holographic description of the multiverse is its potential relevance for the measure problem. 
We have argued that a Wilsonian UV cutoff in the boundary theory corresponds to a scale 
factor cutoff on super-horizon scales in the bulk. The present work has also revealed some 
subtleties in establishing the duality, particularly in the case of bubbles which accelerate 
faster than the Hubble rate of the parent vacuum. Investigation of these issues may lead to 
a more detailed understanding of the boundary theory and its holographic relation to the 
bulk. 
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Appendix: Special Conformal Transformations and their relation to boosts. 

In the flat chart of a dS space of unit radius, with metric given by 

ds^ = ri-\-dri^ + dx^), (-00 < 77 < 0), (81) 
let us onsider the coordinate transformation 

where = {1], x) and = (0, b). Squares such as x"^ and x'^, as well as the scalar products 
below, are constructed with the Minkowski metric, x -y = —x^y^ + x ■ y. We first note that 
since 6° = 0, we have 

4 = -• (83) 

x'^ r]' ^ ' 

Differentiating both sides of ( l82l) and squaring them, we find x'~^dx''^ = x~'^dx'^. It follows 
that (|82l) is an isometry of dS 

7]'-\-d7]'^ + dx'') = r]-^{-drf + dx^). (84) 
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At the future boundary rj = t]' = 0, the coordinate transformation reduces to the Special 
Conformal Transformation (SCT): 



-r. = -.-b'- (85) 



Since fl82l) is an isometry of dS, it must correspond to boosts and rotations in the em- 
bedding Minkowski space with coordinates (T, F, X), where de Sitter space is given by the 
hypersurface X^ + — = 1. Let us see this in exphcit form. Introducing the null 
coordinates U = T — Y and V = T + Y , the embedding coordinates are given by 

U = -a;Vr/, (86) 
V = -Ih, (87) 
X = -x/r/. (88) 



Using (!83l) we have 



^ = ^ = l-26-a; + 6V. (89) 
V x'^ 



where in the last step we have used the square of Eq. (|82|) . After some simple algebra, we 
find 

U' = U, (90) 

V = 1/ - 2b • X + fe^f/, (91) 

X' = X - Uh. (92) 

Hence, in the embedding space, the transformation can be seen as a boost and a rotation, 
which mixes the spatial coordinates X with Y and T. 
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